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Abstract 

We study the effects of sequential fourth quark generation in rare B — > K(K*)t + 1~ 
decays induced by b — > s transition and in B® — B® mixing. Using the experimental 
values on branching ratios of the B — > X s ^f and B — > K(K*)£ + £~ decays, the allowed 
regions for (V^V^I and iVJ'bV^J are determined as a function of the t' quark mass. 
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1 Introduction 



Despite Standard Model (SM) describes successfully all low energy experiments, it is an 
incomplete theory. This theory contains many unsolved and fundamental problems, such as 
the origin of CP violation, mass spectrum and the number of generations. Recent observa- 
tion of neutrino oscillations [1] indicated that neutrino sector of SM must be enlarged. One 
of the most straightforward and economical extension of the SM is adding the fourth gener- 
ation to the fermionic sector, similar to the three-generation case. The extra generation can 
contribute to the electroweak radiative corrections. Latest studies in the electroweak sector 
allow the existence of a fourth generation with heavy Dirac neutrino [2, 3]. Moreover, two 
or three extra generations with relatively "light" neutrinos, with a mass about 50 GeV are 
also allowed [3]. Flavor-changing neutral current (FCNC) transitions provide potentially 
the most sensitive and stringiest test for the SM at loop level, since they are forbidden in the 
SM at tree level. At the same time these transitions are very sensitive to the new physics 
beyond the SM via contributions of the new particles to the loop diagrams. It should be 
stressed that if newly proposed particles are heavy and if they cannot be produced directly 
in the accelerators, their influence through the loop diagrams can be a unique possibility 
for establishing new physics beyond the SM. The effects of the fourth generation to rare 
decays have been studied in many works [4]-[8]. 

Although theoretically FCNC processes are highly suppressed in the SM, very excit- 
ing results are obtained on the experimental side. The first measurement of the FCNC 
processes through b — > S7 were reported by CLEO [9]. Recently, B — > K£ + £~ decay is 
observed at B factories at SLAC and KEK [10]-[12]. BaBar Collaboration also reported 
their preliminary results about the observation of B — > K*£ + £~ decay with branching ratio 
B(B -> K*£ + £~) = (l.68lg;i| ±0.28) x 10- 6 and 90% C.L. B(B -> K*£ + £~) < 3 x lO" 6 
[12]. In this paper, we study the contributions of the fourth generation to the processes 
induced by the b — > s transitions and use the experimental results of the branching ratio 
for the b — > S7 [13], B — > K£ + £~ and B — > K*£ + £~ decays and try to determine the 
constraints on the extended Cabibbo-Kobayashi-Maskawa matrix (CKM) elements (V^V^I 
and \V t , b V; s \. 

The paper is organized as follows. In section 2, we present the basic theoretical ex- 
pressions for the differential widths of the B — > K£ + £~ and B — > K*£ + £~ decays, and for 
the mass difference Am Bs , with sequential up-like quark in the fourth generation model. 
Section 3 is devoted to the numerical analysis and the conclusion. 

2 Theoretical results 

In this section we present the necessary theoretical formulae for the B — > X s 7, B — > K£ + £~ , 
B — > K*£ + £~ decays and for the mass difference in the B® — B® system in presence of the 
fourth generation. All these processes, except mixing in the B® — B® system, are induced by 
the b — > s transition. At quark level this transition is described by the effective Hamiltonian 
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where the full set of operators in Oi(fi) and the corresponding expressions for the Wilson 
coefficients in the SM3 (here and in all further discussions SM3 and SM4 will denote SM 
with three and four generations, respectively) are given in [14, 15]. As is well known, the 
fourth generation introduces the first three generations into SM copiously, and hence it is 
clear that it changes only values of the Wilson coefficients C 7 (/i), C 9 (/x) and Cio(/x) with 
the help of the running fourth generation up quark t' at loop level, and do not introduce 
any new operator structure, i.e., 
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where V t 'b and V t > s are the elements of the 4x4 Cabibbo-Kobayashi-Maskawa (CKM) 
matrix. The explicit forms of the C\ can easily be obtained from the SM results by simply 
substituting m t — > m t >. Neglecting the s quark mass, the effective Hamiltonian leads to the 
following matrix element for the b — > s£ + £~ transition 



M 



Ga 



v tb v* 



cf ^(i - 75)6^7/ + c\p^(i - i 5 )be w e 



- 2Cf^^^(l + 75)^7/ 
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where q 2 = (pi +P2) 2 and pi and P2 are the four-momenta of the final leptons. We observe 
from Eq. (3) that in order to calculate the matrix element for the B — » K*(K)£ + £~ decay, 
the matrix elements of the quark operators in Eq. (3) need to be sandwiched between initial 
and final (K or K*) meson states, which results in a form that is parametrized in terms of 
the form factors 



(K*(p K *,e)\s lt ,(l- l5 )b\B(p B )) = 

-ie*(m B + rriK^A^q 2 ) + i(p B + p K *)n(£*q) 



A2{q 2 ) 
rriB + m,K* 



[As(q 2 ) - A (q 2 ] 



m B + m,K* 
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where e is the polarization vector of K* meson and q = Pb ~Pk* is the momentum transfer. 
Using the equation of motion, the form factor A 3 (q 2 ) can be written in terms of A\(q 2 ) and 
A 2 (q 2 ) as follows 



As(q 2 ) = 



(tub + m K *) 
2rriK* 



A 1 (q 2 )- { "" ! * 1 



2m 



K* 



-A 2 (q 2 ) . 
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In order to ensure that there exists no kinematical singularity we assume that A 3 (q 2 = 0) 
A (q 2 = 0). 
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The corresponding form factors are defined through the matrix elements for the B — > K 
transition as 



(K(p K ) \s^b\ B(p B )) = f + 



m 2 B — m\ 



+ /o 



m%-m K 
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Finetness of Eq. (6) is guaranteed by demanding /+(0) = /o(0). 

The semileptonic form factors for the K* and K mesons resulting from the dipole 
operator sia flv q v {\ + 75)6 are defined as 

(K*(p K .,e) \sia^q»(l + l5 )b\ B(p B )} = 

2e fluXa e*»p x K *q°T l (q 2 ) + 1 [e*(m| - m 2 K .) - (p B + p K *)»(e*q)] T 2 (q 2 ) 
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(K(p K ) \sia^ v q u {l + 75 )6| B(p B )) = - 
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(Pb + Pk)^ 2 ~ - m 2 K ) . (8) 



ffi B + m K 

The matrix elements of the B — > K£ + £~ and B — > K*£ + l~ decays can be written as 
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where i — 1 corresponds to if meson and i = 2 corresponds to X* meson, respectively, and 
1 
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where 
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Using the matrix element, for the dilepton invariant mass distribution we get 

dT— _ G 2 a 2 m% ]VtbV;f ^ J A (3 _ y2) (|Af + |cf) + s(l _ y2){2 + 2r _ s) ]c , f 



ds 
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where v 2 = 1 — Amj/q 2 and A(a, 6, c) = a 2 + b 2 + c 2 — 2ab — 2ac — 26c is the usual triangle 
function. 

In constraining up quark type fourth generation effects, we will also consider B° s — B° s 
mixing. The mass difference Am Bs in SM4 can be written as 
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G 2 m 2 v 
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+ 2rj w (V tb V*) (V t , b V t r s ) S(x u a*)} , (17) 
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where x t = mf/m^, Xf = mt'/ m w an d 
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Here rj t = 0.55 is the QCD correction factor. Taking into account the threshold effect from 
b' quark, r\ tt i is estimated to be [7] 
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Note that when m t i lies between 250 GeV and 400 GeV, r) tt > is quite close to rjf numerically, 
hence for simplicity, in further analysis we will set f] U ' = i] t r . 

In order to obtain quantitative results the value of the fourth generation CKM matrix 
element V t ' b V t f s is needed. For this aim we will use the experimentally measured values 
of the branching ratios BiB — > X s ^i) and BiB — > X c eu e ). To eliminate the uncertainty 
coming from b quark mass we consider the ratio 



R = 
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BiB -> X c ev e ) • 
In leading logarithmic approximation this ratio is equal to 
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where m c = m c /m b and the functions f(m c ) and K(m c ) for the b — * civ transition are given 
by [16] 



f(m c ) = 1 _ 8ml + 8m^ — ml — 2Am\ ln(m c ) 
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From Eqs. (21) and (22) we get 
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The model parameters can be constrained from the measured branching ratios of the B 
X sl decay and B(B -> X c eu e ) = 10.4% 



B(B - X sl ) = 



(3.21 ± 0.43 ± 0.27t{uo) x 10- 4 [17] , 
(3.36 ± 0.53 ± 0.42 ± 0.54) x 10~ 4 [18] , 
[ (3.11 ±0.80 ±0.72) x 10" 4 [19] . 



In our numerical analysis, we will use the weighted average value B(B — > X s 7) = 
(3.23 ± 0.42) x 10" 4 [20] for the branching ratio of the B -> X s 7 decay. 

Another constraint to the extended CKM matrix element comes from the unitarity 
condition, i.e., 
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\v tb 
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1 
1 
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Since charged-current tree-level decays are well measured experimentally they are not 
affected by new physics at leading order. Therefore, for the parameters \V US \, \V CS \, \V cb \ 
and [Ktfe/Vcbl we will make use of Particle Data Group (PDG) constraints [13]. Using the 
weighted average for BiB — > X s ^) and PDG constraint 0.38 < |V c b| < 0.044, from Eqs. 
(24) and (25) we get the constraints 



0.011 < 



C? M V tb V* + Cr w V t , b V t 1 s \ < 0.015 , 
0.03753 < \V tb V* + V t , b V t 1 s \ < 0.043976 , 
< \V ts \ 2 + \V t , s \ 2 < 0.00492 , 
0.998 < |\4| 2 ± |^f,| 2 < 0.9985 . 



(26) 
(27) 
(28) 
(29) 



3 Numerical analysis 

In this section we will study the constraints to | V^V^J coming from the measured branching 
ratios of the B — > K£ + £~ and B — > K*l + l~ decays and B° s — B° s mixing, as well as using 
the results in Eqs. (26)-(29). The main input parameters involved in calculation of the 
branching ratios of the B — > K£ + £~ and B — + K*£ + £~ decays are the form factors, whose 
values we take from light cone QCD sum rule [21]-[23], where the form factors are expressed 
in terms of three parameters as 



1 — CL F S + OpS 

where the values of parameters -F(O), a F and b F for the B — > K and B — > K* decay are 
listed in Table 1. 

The values of the other input parameters which we use in our numerical calculations 
are: m b = 4.8 GeV, m c = 1.35 GeV, m Bs = 5.369 GeV, t Bs = 1.64 x 10" 12 s and 
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yB^K* 
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1 55 


575 


T B^K* 


0.379 


1.59 


0.615 


rpB^K* 

1 2 


0.379 


0.49 


-0.241 


T B^K* 


0.261 


1.20 


0.098 



Table 1: B meson decay form factors in a three-parameter fit, where the radiative correc- 
tions to the leading twist contribution and SU(3) breaking effects are taken into account 
(see [22, 23]). 

B B J Ba = (0-26 GeV) 2 . The experimental lower bound of the mass difference is Am Bs > 
14.9 ps^ 1 . For the values of the Wilson coefficients Cf M , Cg M and Cf M we have used their 
next-to-leading logarithmic result: Cf M = -0.308, C 9 5M = 4.154 and Cf M = -4.261. 
It should be noted that the decays B — > K£ + £~ and B — > K*£ + £~ receive long distance 
contribution coming from cc intermediate states. In the present work we neglect such long 
distance effects. The strong dependence on m t > (see for example [4]) makes the electroweak 
penguins a good place looking for the existence of fourth generation. Contributions of fourth 
generation to B — ■> K(K*)£ + £~ decays have already been studied (see second references in 
[4] and [6]). The present investigation differs from the above-mentioned works in two 
aspects: 

• we use the experimentally measured results on branching ratio, 

• we consider V t bV t * and V t >bV t * s as two independent complex parameters, which were 
taken to be real in [4] and [6]. 

The complex parameters V t bV t * and V t >bV t * s are constrained by the unitarity conditions 
(see Eqs. (27)-(29)), the measured branching ratios B — > X s 7 (see Eq. (26)) and B — > 
K(K*)£ + £~ decays (see Eqs. (15), (16)), which depend on m t i. For each value of m t i there 
exists an allowed region in the |Vt&V^*|-|Vt'&V^ s | plane. 

Since there exists no analytical solution of Eqs. (26)-(29), we will solve these equations 
numerically assuming that V t bV t * s and V t 'bV t J s are complex. For sufficiently large number of 
randomly chosen complex parameters Vt'bV t * s and V t bV t *, the selected values would range 
over the whole solution space. In Figs. (l)-(3) we present the allowed region for VtbV t * s and 
V t 'bV t 7 s at m t i = 200 GeV, m t > = 300 GeV and m t i = 400 GeV, respectively. In obtaining 
this solution region we have used Eqs. (26)-(29) and Eq. (15). From these figures we see 
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that, for iVj/tV^J = 0, |Vt&VJ*| takes on values that is close to the SM prediction and is 
mainly distributed around ~ 0.04. Moreover, when |V^ftV^,| increases, the allowed region of 
|VtbV^*| becomes wider with the center being fixed around 0.04, and the values of |Vt'bV^ s | 
are mainly distributed around 0.01. With increasing values of m t ', the allowed region for 
IK'&V^J becomes narrower. This behavior can be explained as follows. Wilson coefficients 
are strongly dependent on m t > and in order to remain in the experimentally allowed region, 
the element |V^VJ*J of CKM matrix must decrease, since branching ratio contains factors 
like |CjV#,V^*|. Similar behavior is observed for the B — > K*£ + £~ decay (see Figs. (4)-(6)). 

In Figs. (7) and (8) we present the dependence of Am Bs on (V^V^I and IV^V^J, taking 
into account the lower experimental bound for ArriB s , at two different values of m t t. It 
follows from both figures that the main distribution is in the range 0.36 < | Vjj, V^.* | < 0.044 
and < |Vt'&V^*J < 0.01. With increasing values of m t >, obviously, Atob s also increases. 

Finally we would like to note that restrictions to |Vf/dV^ 6 | matrix element can be obtained 
by an investigation of the rare decays induced through b — > d transition. Further analysis 
of the decays induced by b — > s(d) transition is more promising in studying new sources for 
CP violation, since 4x4 CKM matrix predicts the existence of three CP violating phases. 
We will discuss this issue elsewhere in future. 

In conclusion, we have studied the effect of the fourth generation quark to the rare 
decays induced by FCNC b — > s transition. Using the experimental result for the branching 
ratios of the B — > X s 7, B — > K(K*)l + l~ decays and the unitarity condition for the 4x4 
CKM matrix, we have determined the allowed parameter space for (V^V^I and |VWV^*J m 
their dependence on m t '. 
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Figure captions 



Fig. (1) Three-dimensional plot of the branching ratio for the B — > K£ + £ decay, with 
respect to the allowed parameter space of V t bV t * s and V t 'bV t * s , at m t > = 200 GeV. 
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Fig. (7) Three-dimensional plot of the mass difference Am# s of the B® — B® system, 
with respect to the allowed parameter space of V t bV t * and V t 'bV t * s , at m t > = 300 GeV. 

Fig. (8) The same as in Fig. (7), but at m v = 700 GeV. 
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